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1.0  INTRODUCTION 

Renewed  interest  in  departure  characteristics  of  hi^  performance  military  aircraft 
over  the  past  decade  has  resulted  in  the  development  of  a  captive  aircraft  departure 
system  (CADS)  testing  technique  at  the  Arnold  Engineering  Development  Center  (AEDC). 
The  CADS  testing  technique  (Ref.  1)  is  capable, of  generating  aircraft  maneuvers  up  to 
and  including  departure  in  the  wind  tunnel  utilizing  the  tunnel  as  an  analog  forcing 
function.  Aircraft  static  external  forces  and  moments  required  in  the  computer  simulation 
are  measured  during  the  course  of  the  maneuver  with  all  dynamic  components  being 
input  parameters.  In  the  past,  when  aircraft  operating  envelopes  were  in  the  linear 
aerodynamic  (low  angle-of-attack)  flight  regime  these  dynamic  components  consisted 
primarily  of  direct  and  cross-damping  derivatives.  At  the  high  angle-of-attack,  severe 
maneuvering  conditions  of  modem  figliter  aircraft,  there  is  reason  to  believe  that 
additional  dynamic  derivatives  (cross  coupling)  may  be  needed  to  insure  accurate  motion 
simulation.  This  assumption  is  based  on  work  reported  by  Orlik-Ruckcmann  and  others  in 
Ref.  2.  Experimental  data  in  Ref.  2  indicate  that  at  high  angles  of  attack  (maneuvering 
flight  condition)  both  cross  and  cross-coupling  dynamic  derivatives  acquire  values  orders 
of  magnitude  larger  than  those  at  the  lower  angles  (cruise  flight  condition).  The  subject 
report  investigates  the  effect  of  these  large  values  in  cross  and  cross-coupling  derivatives 
on  aircraft  motion  and  thereby  identifies  the  important  dynamic  derivatives  to  be 
included  in  the  CADS  and  similar  motion  simulation  studies. 

Aircraft  motion  consists  of  combinations  of  the  basic  longitudinal  and  lateral 
directional  characteristic  modes  of  motion.  The  sensitivity  of  each  of  these  characteristic 
modes  to '  variations  in  the  dynamic  derivatives  is  an  indicator  of  the  derivatives  wliich 
should  be  included  in  motion  simulation  studies.  A  5  deg  of  freedom  (5-DOF)  linearized 
stability  program  is  utilized  in  this  report  for  identifying  mode  changes  with  variations  in 
derivatives.  The  linearized  prop^m  identifies  the  basic  aircraft  modes  with  the  exception 
of  the  longitudinal  phugoid.  Simplifying  assumptions  used  in  linearizing  the  program 
required  that  the  aircraft  velocity  remain  constant,  thereby  eliminating  the  phugoid 
mode.  To  simulate  maneuvering  fliglit,  the  linearized  analysis  is  performed  in  both  level 
and  turning  flight.  Changes  in  the  motion  modes  with  derivative  variations  are  depicted  in 
root  locus  diagrams.  Three  typical  fighter  aircraft  configurations  are  investigated  in  the 
analysis,  the  F-4,  A-7,  and  F-5  aircraft.  These  aircraft  were  chosen  because  each  is 
characterized  by  distinct  flying  qualities  which  represent  areas  of  interest  in  flight 
dynamics  of  modern  high  performance  military  aircraft.  In  addition  to  the  above  cross 
and  cross-coupling  derivative  investigation,  the  importance  of  the  direct  acceleration 
derivative,  C„p,  in  high  angle-of-attack  motion  simulation  is  reviewed. 
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2.0  BACKGROUND 

Aircraft  dynamic  damping  derivatives  are  normally  categorized  into  three  distinct 
groups:  direct,  cross,  and  cross-coupling.  The  primary  derivatives  associated  with  each  of 
these  groups  are  shown  in  Fig.  1.  At  low  angles  of  attack,  which  are  characteristic  of 
aircraft  cruise  flight  conditions,  the  dynamic  cross  and  cross-coupling  derivatives  can  be 
predicted  analytically  (Refs.  3  through  6).  Calculations  using  the  prediction  techniques 
normally  produce  quite  accurate  results  in  both  longitudinal  and  lateral  directional 
derivatives  (Ref.  7).  These  derivatives  normally  exhibit  small  variations  with  changes  in 
flight  conditions  at  low  angles  of  attack,  thereby  producing  a  small  constant  effect  on  the 
aircraft  flight  characteristics.  As  a  result,  derivatives  are  often  input  to  motion  simulation 
programs  as  constant  values  with  good  success.  Both  calculated  and  measured  values  of 
the  cross-coupling  derivatives  in  this  flight  regime  are  small  and  are  normally  omitted 
from  the  motion  simulation. 

On  the  other  hand,  the  flight  envelope  of  the  current  generation  of  fighter  aircraft 
necessitates  improved  prediction  and  motion  simulation  techniques  because  of  the 
aircraft's  greatly  expanded  operational  envelope.  The  high  angle  of  attack  and  high  speeds 
at  which  these  aircraft  operate  produce  nonlinear  aerodynamic  flow  phenomena  sudh  as 
separated  flows,  vortex  shedding,  etc.  Not  only  these  phenomena,  but  also  certain 
asymmetric  flows  such  as  separation  on  the  leeward  wing  or  vortex  shedding  from  the 
aircraft  forebody,  which  can  occur  at  high  angles  of  sideslip,  all  have  a  significant 
influence  on  the  aircraft  dynamic  derivatives.  The  cross-coupling  derivatives  previously 
considered  to  have  values  near  zero  assume  values  which  are  10  to  15  percent  of  the 
cross  derivatives  on  current  fighter  configurations  (Ref.  2).  The  cross  and  cross-coupling 
derivatives  not  only  increase  in  magnitude  but  often  experience  sign  changes  over  very 
small  angle-of-attack  ranges.  Because  of  these  effects,  all  dynamic  derivatives  may  be 
required  inputs  in  motion  simulation  programs  as  variables,  normally  functions  of  angle 
of  attack  and  angle  of  sideslip. 

All  derivatives  listed  in  Fig.  1  may  not  be  required  for  achieving  accurate  motion 
simulation.  The  direct  damping  derivatives  are  always  considered  of  primary  importance 
and  are  included  in  motion  simulation  studies  at  low  and  high  angles  of  attack.  This 
report  addresses  the  importance  of  the  remaining  cross  and  cross-coupling  derivatives  at 
high  angles  of  attack. 


3.0  AIRCRAFT  CONFIGURATIONS 

Three  aircraft  configurations  were  chosen  for  this  dynamic  sensitivity  study.  Each 
aircraft  was  selcted  on  the  basis  of  inherent  high  angle-of-attack  lateral-directional  flight 
characteristics.  The  intent  was  to  select  aircraft  which  would  exhibit  a  range  of 
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lateral-directional  stability  characteristics  typical  of  current  high  performance  military 
aircraft.  On  the  basis  of  these  criteria,  configurations  similar  to  the  F^,  A-7,  and  F-5 
aircraft  were  selected  for  the  analysis. 

The  F-4  configuration,  a  twin-jet,  swept-wing  fighter  aircraft,  is  shown  in  Fig.  2. 
This  configuration  possesses  a  "wing  rock"  or  dutch  roll  oscillation  at  high  angles  of 
attack  (Refs.  8  and  9).  As  angle  of  attack  increases  up  to  approximately  23  deg,  the 
dutch  roll  oscillation  becomes  unstable,  resulting  in  aircraft  directional  divergence. 

The  A-7  configuration,  a  single-engine,  swept-wing  aircraft,  is  shown  in  Fig.  3.  This 
configuration  exhibits  some  directional  instability  at  angles  of  attack  above  20  deg  while 
still  possessing  lateral  stability  (Refs.  9  and  10).  The  A-7  departure  is  characterized  by  an 
abrupt  "nose  slice". 

The  F-5  configuration,  a  twin-engine,  tapered-wing  aircraft,  is  shown  in  Fig.  4.  This 
configuration  exhibits  excellent  lateral-directional  dynamic  stability  characteristics  at 
angles  of  attack  well  above  that  required  for  wing  stall  (Ref.  II).  The  wing-rock  and 
nose-slice  phenomena  experienced  by  the  F^  and  A-7  configurations  are  not  experienced 
in  the  F-5  configuration. 


4.0  METHOD  OF  ANALYSIS 

The  sensitivity  of  aircraft  motion  to  variations  in  aerodynamic  and  physical 
characteristics  is  reflected  in  its  lon^tudinal  and  lateral-directional  stability  modes,  the 
two  longitudinal  modes  being  short  period  and.  phugoid,  and  the  three  lateral-directional 
modes  being  roll,  spiral,  and  dutch  roll.  These  characteristic  modes  are  normally 
investigated  independently  through  linearized  3  deg  of  freedom  (3-DOF)  stability 
programs.  Such  programs  provide  the  damping  and  frequency  characteristics  of  the  modes 
with  the  aircraft  in  level  flight,  with  zero  angular  rates  of  motion.  However,  the  flight 
regimes  addressed  in  this  report  consider  both  level  and  turning  flight  and  therefore 
cannot  be  correctly  investigated  with  only  the  3-DOF  program.  In  addition,  the 
longitudinal  and  lateral-directional  flight  modes  are  separate  in  the  3-DOF  analysis; 
therefore,  the  cross-coupling  dynamic  derivatives  are  nonexistent  in  the  separated  modes. 

A  5-DOF  linearized  stability  program  has  been  developed  for  conducting  the 
dynamic  sensitivity  study.  In  this  program,  the  reference  flight  condition  is  generalized  to 
include  the  effect  of  non-zero  bank  angle,  pitch  rate,  and  yaw  rate.  The  result  is  a  weak 
coupling  between  longitudinal  and  lateral-directional  dynamics  providing  the  mechanism 
for  investigating  the  importance  of  the  cross-coupling  dynamic  derivatives  in  turning 
flight.  The  5-DOF  system  was  obtained  by  eliminating  the  X-body  axis  translation 
equation  from  the  six  Euler  equations  of  motion.  Omitting  this  equation  requires  that  the 
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aircraft  remain  at  a  constant  velocity  along  the  X  axis,  thereby  alleviating  the  longitudinal 
phugoid  mode  of  motion.  The  Uneraized  equations  of  motion  derived  with  these  criteria 
are  listed  below.  The  linearization  process  used  is  presented  in  Refs.  3  and  4. 


-  N^/3  -  =  0 


(5) 


0  —  cos  0^4  +  sin  (fj =  0 


(6) 


0  ~  P  - 


where 


"i  - 


-  I 


yy 


1  -  I 


yy 


-  I. 


(7) 


Equations  (1)  and  (2)  are  the  force  equations;  (3),  (4),  and  (5)  are  the  moment 
equations;  and  (6)  and  (7)  the  kinematic  equations  which  the  motion  must  satisfy.  Initial 
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pitch  attitude  and  initial  roll  rate  are  assumed  to  be  small  in  the  derivation.  The 
equations  are  referenced  to  the  Eulerian  body  axis  system  with  the  axes'  origin  at  the 
aircraft  center  of  gravity.  The  aircraft  XZ  plane  is  the  plane  of  symmetry. 

Initial  inputs  q^,  and  r,,,  required  for  solving  the  equations,  are  evaluated 

from  the  conditions  for  steady,  turning  flight. 


n  '^\n 

^0 

~  V  ® 

so 

=  sin 

(9) 

=  ijj^  cos 

(10) 

By  applying  the  Laplace  operator,  one  obtains  the  7x7  characteristic  determinant  shown 
in  Appendix  A.  Values  ei  through  ei4  of  the  determinant  correspond  to  terms  in 
parentheses  in  the  respective  equations  (1)  through  (7).  Solving  this  determinant  yields 
the  characteristics  equation  below,  with  coefficients  also  defined  in  Appendix  A. 


f.  ,|A'  -  +  tAC2  +  ABl  ^  +  (AB2  •  AC3  +  <A„ADi  ^ 

f  +  (AB3  +  AC4  1-  dy„AD2  -  d'„AE2  +  v>„AC2)A^ 

f  (AB-t  ^  AC5  t  H  -  ^i.“AC3)A“ 

+  IAH.3  +  w^ADt  -  (/lyAKq  i  i.!'|“ACriA 

4  AD.:,  +  ACS  +  w;AC,->!  .0  (11) 


This  algebraic  expression  for  the  characteristic  determinant  (Appendix  A)  was  formulated 
for  use  m  future  work  of  generating  aircraft  stability  boundaries  by  utilizing  the  Taylor 
series  expansion  technique  of  Ref.  12 

With  tlie  seventh-order  characteristic  equation,  three  new  eigenvalues  appear  in 
addition  to  those  which  normally  occur  with  a  3-DOF  lateral-directional  system.  Two  of 
the  eigenvalues  represent  the  longitudinal  short  period  mode,  and  the  third  eigenvalue 
represents  a  coupling  with  the  lateral-directional  spiral  mode  to  form  an  oscillatory  spiral. 
The  remaining  lateral-directional  modes  are  the  conventional  dutch  roll  and  the  roll. 
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5.0  RESULTS  AND  DISCUSSION 


5.1  GENERAL 

The  aerodynamic  and  physical  characteristics  for  each  of  the  aircraft  configurations 
evaluated  are  presented  in  Tables  1,  2,  and  3. 

The  F-4  characteristics  are  given  in  Table  1.  The  mass,  inertia,  and  geometric 
parameters  shown  represent  a  clean  configuration  as  outlined  in  Refs.  7,  13,  and  14. 
The  aerodynamic  static  and  dynamic  stability  derivatives  were  acquired  from  Refs.  7  and 
13.  The  longitudinal  and  lateral-directional  derivatives  were  obtained  over 
angle-pf-attack  and  angle-of-sides!ip  ranges,  respectively,  of  ±5  deg  about  the  aircraft  trim 
angle  of  attack  of  29  deg.  The  trim  angle  of  attack  represents  the  maximum  trim 
capability  of  the  horizontal  tail  (it  =  -25)  based  on  static  force  tests. 

The  A-7  characteristics  are  shown  in  Table  2.  The  mass,  inertia,  and  geometric 
parameters  are  for  a  clean  aircraft  configuration  (Ref.  15).  All  aerodynamic  characteristics 
are  from  unpublished  static  wind  tunnel  data.  Static  derivatives  were  calculated  over 
angle-of-attack  and  angle-of-sideslip  ranges  of  ±5  deg  about  a  nominal  trim  angle  of 
attack  of  26.5  deg.  The  26.5-deg  attack  angle  corresponds  to  a  maximum  horizontal 
stabilizer  setting  of  -25  deg  (it  =  -25). 

The  F-5  aerodynamic  characteristics  are  shown  in  Table  3.  Mass,  inertia,  and 
geometric  characteristics  are  representative  of  those  for  the  clean  F-5E  aircraft.  All 
aerodynamic  characteristics  were  acquired  from  wind  tunnel  data  presented  in  Ref.  11. 
Both  static  and  dynamic  derivatives  were  evaluated  about  an  angle’  of  attack  of  28  deg 
with  static  derivative  calculations  being  over  a  ±5-deg  range.  The  28-deg  angle  of  attack 
corresponds  to  an  unrealistic  horizontal  stabilizer  setting  of  -25  deg  (q  =  -25).  Maximum 
horizontal  stabilizer  deflection  on  the  F-5  aircraft  is  -17  deg.  The  increased  deflection  was 
used  to  produce  higher  angles  of  attack  representative  of  flight  at  load  factors  greater 
than  Ig. 


The  dynamic  stability  analysis  on  all  three  aircraft  was  conducted  in  1-g  (level),  3-g 
(turning),  and  6-g  (turning)  flight  at  an  altitude  of  30,000  ft.  The  turning  flight 
conditions  simulate  the  aircraft  in  (he  maneuvering  flight  regime.  To  isolate  the  effects  of 
the  longitudinal  and  lateral-directional  coupling,  which  occur  in  turning  flight,  on  the 
dynamic  sensitivity  study,  it  was  deemed'  necessary  to  maintain  identical  aerodynamic 
stability  derivatives  at  all  three  flight  conditions.  This  was  accomplished  by  maintaining  a 
constant  angle  of  attack  at  all  flight  conditions,  The  aircraft  airspeed  was  adjusted  to 
achieve  the  desired  load  factor  and  maintain  a  constant  altitude.  It  was  assumed  that  the 
aerodynamic  coefficients  do  not  change  significantly  over  the  speed  ranges  encountered. 
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A  summary  of  the  flight  conditions  for  simulating  the  aircraft/flight  combinations  utilized 
is  given  in  Table  4. 

5.2  DYNAMIC  CROSS  DERIVATIVE  AND  C„^  VARIATIONS 

Fighter  aircraft  motion  sensitivity  to  variations  in  the  dynamic  cross  derivatives 
Cnp,  and  is  shown  in  Figs.  5,  6,  and  7,  respectively.  Each  figure  represents  changes  in 
the  F-4,  A'7,  and  F-5  aircraft  longitudinal  and  lateral-directional  modes  of  flight  with 
variations  in  the  dynamic  cross  derivatives.  As  noted  previously,  the  longitudinal  phugoid 
motion  is  not  included.  Each  cross  derivative  is  varied  over  a  range  of  +5  per  radian,  a 
variation  which,  in  most  cases,  is  an  order  of  magnitude  greater  than  that  nominally 
e.xperienced  by  fighter-type  aircraft.  As  the  derivatives  are.  varied,  the  resultant 
eigenvalues  of  the  characteristic  equation  are  plotted  in  the  root  locus  format  shown  (real 
and  imaginary  axes  representing  damping  and  frequency,  respectively).  The  time  required 
for  each  mode  to  damp  to  one  half  or  diverge  to  twice  amplitude  is  shown  horizontally 
below  each  plot  with  the  motion  period  being  shown  vertically. 

The  longitudinal  and  lateral-directional  stability  modes  of  each  aircraft  are  presented 
for  1-g  level  and  3-g  and  6-g  turning  flight.  In  all  cases,  the  motion  sensitivity  to  Cg^., 
Cjip,  and  variations  increased  with  the  increasing  load  factor  or,  in  this  case,  the 
steady-state  turn  rate. 

When  transitioning  from  level  to  turning  flight  the  aircraft  spiral  mode  becomes 
oscillatory.  The  additional  eigenvalue  causing  the  oscillatory  spiral  results  from  the 
non-zero  turning  rate  (ij/o)  term  in  the  characteristic  equation  |Eq.  (11)].  With  zero 
Eq.  (11)  is  sixth  order  and  the  spiral  mode  uncoupled.  The  last  two  coefficients  of  Eq. 
(11)  represent  the  spiral  flight  mode.  The  total  absence  of  the  F-5  spiral  mode  in  level 
flight  results  from  a  combination  of  zero  i/Zq,  Cn^,  and  Cgp  terms.  These  zero  terms  cause 
the  last  two  coefficients  of  Eq.  (11)  to  vanish,  resulting  in  the  fifth-order  characteristic 
equation.  Zero  values  of  Cg^  and  Cn^  correspond  to  measured  damping  derivatives  (Ref. 
1 1)  occurring  at  the  trim  flight  condition  selected  (Table  3). 

Yaw  Rate  Parameter  ■  Both  F-4  and  A-7  roll  and  dutch  roll  modes  of  flight  (Fig.  5) 
show  sensitivity  to  variations  in  the  rolling  moment  cuased  by  yaw  rate  parameter  Cg^. 
The  dutch  roll  mode  of  the  F-4  goes  from  a  divergent  oscillation  with  positive  values  of 
Cgj.  to  an  aperiodic  and  finally  roll-coupled  mode  at  negative  values.  The  roll  and  dutch 
roll  modes  of  the  A-7  may  be  highly  divergent  or  stable,  depending  on  the  magnitude  and 
sign  of  Cg^.  The  F-5  stability  modes  show  virtually  no  sensitivity  to  Cg^  variations. 

Roll  Rate  Parameter  -  Variations  of  aircraft  yawing  moment  caused  by  roll  rate 
parameter  Cpp  affect  damping  of  the  F-4,  A-7,  and  F-5  dutch  roll  and  roll  modes  of 
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flight  significantly,  as  seen  in  Fig.  6.  The  longitudinal  short  period  mode  of  the  F-4  also 
shows  a  reduction  in  damping  approaching  divergence  as  C„p  becomes  more  negative.  For 
most  aircraft,  the  more  negative  C„p  is,  the  greater  is  the  reduction  in  dutch  roll 
damping.  This  trend  exists  in  the  A-7  and  F-5  aircraft,  but  not  in  the  F-4,  where  the 
dutch  roll  mode  actually  becomes  stable  for  values  of  C„p  less  than  approximately  -2. 

Dynamic  Stability  Derivatives  -  When  wind  tunnel  testing  is  employed  for  generating 
dynamic  stability  derivatives,  the  and  C„g  terms  are  normally  coupled  with  the 
and  Cnj  terms,  respectively,  in  the  data  acquisition.  As  a  result,  motion  simulation 
programs  often  combine  these  terms  as  (Cg^  +  Cgp  and  (C„j  +  C„^)  and  operate  on  them 
with  the  angular  rate  term  r.  This  technique  assumes  the  predominant  term  in  the 
expressions  to  be  Cg^  and  Cn^.  Recent  experimental  data  (Ref.  16)  have  shown  fallacies 
in  this  assumption  during  operation  at  high  angles  of  attack.  At  high  angles  of  attack,  the 
Cg^  and  Cn^  terms  become  predominant,  with  the  Cg^  and  Cn^  terms  approaching  zero. 
The  effects  that  large  positive  and  negative  values  of  Cg^  and  Cn^  may  have  on  fighter 
aircraft  stability  are  shown  in  Figs.  7  and  8,  respectively.  Damping  of  both  the  F-4  roll 
and  the  aperiodic  dutch  roll  roots  of  Fig,  7  is  sensitive  to  positive  and  negative  variations 
in  Cg^.  The  dutch  roll  mode  of  tha  A-7  aircraft  experiences  changes  in  both  frequency 
and  damping  with  Cg^  variations.  As  Cg^  appraoches  2,  the  mode  frequency  approaches 
zero  and  motion  becomes  aperiodic.  The  F-5  dutch  roll  motion  experiences  large  changes 
in  damping  and  frequency  with  Cg^  variations.  The  remaining  modes  of  the  F-5  are  all 
insensitive  to  Cg^  variations. 

Root  movements  identical  to  those  experienced  in  Fig.  7  are  seen  in  the  Cn^ 
variation  of  Fig.  8.  The  sensitivity  of  the  F-4,  A-7,  and  F-5  lateral-directional  stability 
modes  to  variations  of  Cn^  is  not  as  great  as  that  seen  in  the  Cg^  variation  of  Fig.  7  but 
is  still  of  a  magnitude  to  be  considered  in  motion  simulation  work. 

5.3  DYNAMIC  CROSS-COUPLING  DERIVATIVE  VARIATIONS 

The  dynamic  cross-coupling  derivatives  investigated  in  this  analysis  are  Cmp,  C^,, 
Cg;,  Cn^,,  and  Cg^^.  Sensitivity  of  the  F-4,  A-7,  and  F-5  longitudinal  and 

lateral-directional  stability  modes  to  variations  of  these  derivatives  over  ranges  of  ±1  per 
radian  are  presented  in  Figs.  9,  10,  and  11.  The  range  over  which  the  derivatives  are 
varied  was  selected  from  experimental  wind  tunnel  data  obtained  with  an  airplane-type 
configuration  (Ref.  2).  In  all  cases,  the  derivative  range  employed  is  more  than  an  order 
of  magnitude,  and  in  most  cases,  two  orders  of  magnitude,  larger  than  the  range 
exhibited  by  the  experimental  data. 

Cmp.  Cg^,  "  The  group  of  derivatives  shown  in  Fig.  9  all  have  one  thing  in 
common;  they  have  no  effect  on  the  F-4,  A-7,  and  F-5  longitudinal  and  lateral-directional 
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Stability  modes  in  either  level  or  turning  flight.  Based  on  these  results,  the  derivatives 
should  not  be  required  for  acquiring  accurate  motion  simulation  of  fighter  aircraft  in  the 
maneuvering  flight  regime. 

Cn  )  Cc  -  Changes  occurring  in  the  F-4,  A-7,  and  F-5  stability  modes  with  varying 

and  Cgq  are  presented  in  Figs.  10  and  11,  respectively,  These  derivatives  are 
separated  from  the  above  cross-coupling  derivatives  because  they  do  produce  small 
changes  in  some  of  the  aircraft  stability  roots.  The  greatest  root  change  occurs  in  the  F-4 
spiral  and  dutch  roll  modes  with  varying  Cg^.  Because  was  varied  over  a  range  (-1,  1) 
which  is  approximately  two  orders  of  magnitude  greater  than  its  nominal  value  (Ref.  2), 
the  change  in  the  stability  roots  is  considered  insignificant. 

6.0  CONCLUDING  REMARKS 

As  a  result  of  this  analysis,  the  following  observations  and  conclusions  are  offered; 

1.  A  5  deg  of  freedom  linearized  stability  program  provides  a  technique  for 
investigating  lateral-directional  stability  modes  of  an  aircraft  in  level  and 
turning  flight  with  light  longitudinal  coupling. 

2.  When  ascertaining  the  importance  of  a  dynamic  derivative  in  aircraft 
stability  calculations,  one  should  examine  a  number  of  different  aircraft 
configurations;  even  so,  a  general  conclusion  may  be  impossible. 

3.  In  the  5  deg  of  freedom  linearized  analysis,  the  aircraft  stability  modes 
become  more  sensitive  to  derivative  variations  as  the  load  factor  or  turning 
rate  of  the  aircraft  increases. 

4.  The  dynamic  cross  derivatives  Cnp,  Cg^,  and  Cgp  should  be  included  in 
motion  simulation  of  fighter  aircraft  examined  in  the  subject  analysis. 

5.  All  cross-coupling  derivatives  Cg^^,  Cn^,  C^p,  Cm^,  Cg^,  and  C„^  are 
considered  unimportant  in  motion  simulation  of  all  fighter  aircraft 
examined  in  this  analysis. 

6.  The  rate  of  change  of  sideslip  derivatives  Cg^  and  C„^  produces,  in  each 
case,  significant  changes  in  the  lateral-directional  stability  modes  of  the 
aircraft  examined.  Uncoupling  of  the  (Cg^  +  Cgp  and  (Cn^  +  Cnp  terms 
may  be  necessary  for  accurate  aircraft  motion  simulation  during  operation 
at  high  angles  of  attack. 
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Figure  10.  Continued. 
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c.  F-5 

Figure  10.  Concluded. 
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a.  F-4 

Figure  11.  Locus  of  roots  with  variation. 
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Figure  11.  Continued. 
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Table  1.  F-4  Aircraft  Physical  and  Aerodynamic  Characteristics 


Weight 

1,215  slugs 

Ci, 

0,  638 

29,950  slug-ft^ 

Cn|3 

-0. 252 

% 

164,300  slug-ft^ 

-0.031 

169,  538  slug-ft^ 

Cnr 

-0.444 

Ixz 

5,241  slug-ft^ 

Cie, 

0 

S 

538.34 

Ciq 

0 

b 

38.41  ft 

Cmq 

-4.0 

c 

16.04  ft 

Cm 

111^ 

-0.5 

eg 

33  percent  MAC 

Cm- 

0 

a 

29  deg 

Cmr 

0 

Cl 

0.  928 

Cn^ 

0 

-0.201 

cn. 

1.  76 

Cyp 

0.  056 

ca 

0.05 

Cyr 

-0.708 

0 

-0. 017 

Cn^ 

0 

-0.410 

Cnq 

0 
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Table  2.  A-7  Aircraft  Physical  and  Aerodynamic  Characteristics 


Weight 

719  slugs 

0. 

075 

Ix 

15 

,927  slug-fr 

-0. 

123 

V 

64 

.792  slug-ft^ 

o 

H 

0. 

0175 

Iz 

75 

,976  slug-ff^ 

Cn 

“r 

-0. 

23 

^xz 

3, 

885  slug-ft^ 

0 

S 

■37 

5.00 

0 

b 

38 

.73  ft 

Cmq 

-4. 

9 

c 

10 

.84  ft 

Cm 

-1. 

56 

a 

eg 

30 

percent  MAC 

Cm^ 

0 

a 

26 

.  5  deg  ' 

Cm^ 

0 

Cl 

1. 

07 

0 

-0. 

974 

Cn^- 

0. 

871 

0. 

14 

Ca 

-0. 

04 

Cy^ 

0. 

35 

Ci^ 

0 

* 

-0. 

10 

Cn- 

“Of 

0 

Cje 

p 

-0. 

01 

Cn 

"q 

0 
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Table  3.  F-5  Aircraft  Physical  and  Aerodynamic  Characteristics 


Weight 


V 

S 

b 

c 

eg 

Cl 

Cy^ 


Cn 

Cje, 


357  slugs 
3,400  slug-ft2 
37,  BOO  slug-ft^ 
41, 000  slug-ft^ 
-200  slug-ft^ 
181.00  ft^ 

25.80  ft 
7 . 88  ft 

18  percent  MAC 
28  deg 
1.08 
-0.  29 
0.  20 
1.2 
-0. 18 
0 


C., 

-0.05 

0 

0. 175 

Cn^ 

0 

Cje . 

0 

Ci_ 

0 

4 

Cm_ 

-18.4 

-q 

Cm 

lilQ, 

-1.42 

Cm ' 

0 

Cm 

0 

r 

Cn^ 

0 

Cn 

2,98 

a 

ca 

-0.02 

0 

Cpi 

0 

Cnq 

0 
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Table  4.  Flight  Conditions 


deg 

V.. 

ft/sec 

rad/ sec 

rad/ sec 

e> 

rad/sec 

F-4  AIRCRAFT 

1 

0 

78.  3 

419,  2 

0 

0 

0 

3 

70.  5 

238.8 

732.  3 

0. 124 

0. 117 

0.041 

6 

80.4 

469. 6 

1,027.9 

0. 185 

0. 182 

0.031 

A-7  AIRCRAFT 

1 

0 

57.  6 

359.  8 

0 

0 

0 

3 

70.  5 

176.4 

629,  4 

0.  145 

0.137 

0.048 

6 

80.4 

353,0 

890.3 

0.214 

0.211 

0,036 

F-5  aircraft 

1 

0 

58.  8 

363. 0 

0 

0 

0 

3 

70.  5 

173.0 

623.3 

0.  146 

0. 138 

0.049 

6 

80.4 

345.9 

881.3 

0.216 

0.213 

0.036 
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APPENDIX  A 

5-DOF  CHARACTERISTIC  DETERMINANT 


(A  -  Z^)o 

0 

0 

3 

0 

0 

0 

(A  - 

-«iP 

0 

0 

-«2>^ 

-L^a 

—  (ALjg  +  Ljg)  /S 

(A  -  e-)p 

efl'q 

I 

-  (A  -  eg)r 

J(  X 

0 

0 

0 

ejo'P 

(A-M/q 

0 

n 

=  0 

1 

■ 

-(AN^^N^)^  -U^_e,2)p 

*  Z7 

'ja*! 

(A  +  *14^' 

0 

n 

0 

0 

0 

-cos 

sin  i^„r 

\e 

0 

0 

“P 

0 

0 

A<^ 

Coefficients  of  Characteristic  Equation  [Eq. 

(11)] 

ABl  = 

-e^bu  +  egbgj 

AB2  = 

-  Cgbjj)  - 

®2^^3]^\v 

J-  b32  ~  b^,) 

AB3  = 

®5^^v*^12  “  ^13  " 

®3^22^  “  ® 

2*'^w^2  “ 

*^33  *^42^ 

AB4  = 

“  ^14  “ 

®3^23^  “  ® 

2^^w*’33  “ 

^4  +  *^43'^ 

AB5  = 

-  e2(b3^Z^  + 

ACl  = 

1 

n 

+ 

2.^  -  C12 

+  C21 

"  ^31®3  ~  ^41 

AC2  = 

V 

- 

^12^^  V  + 

ZJ  + 

Cj2  +  BjCCgj  -  - 

-  C5,) 

+ 

-  ^32 

-  ^61^  +  ^4l^v  ”  ^42 

AC3  = 

V 

- 

Zw)  - 

Ci4  -  ej(C23  “  C22Z,,,  - 

-  C52) 

+ 

*3^^32^vv 

r 

'-33 

-  0^2)  +  ^42^  V  “  ^43 

AC4  = 

- 

Cu<^'v  - 

z.;  - 

®/’-24  ~  ^23Zvv  ~  ^53^ 

-t- 

“  ^34 

-  C^3)  +  ^43^v  ~  ^44 

ACS  = 

CuZ.Yv 

“ 

®1^^242w 

-  C54) 

+  e3(C34Z^^,  -  Cg^)  +  C 

44^  V 
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ADI  =  (DBBl)  cos  (h  e. 

5 

AD2  =  DBAl  +  e3[(DBB2)  cos  0*,^  _  (DBBl)  Y,,  cos  (^>J  +  DBCl  DBDl 

AD3  =  DBA2  -  e5[(DBB3)  cos  -  (DBB2)Y^  cos  0^]  +  D8C2  +  DBD2 

AD4  =  DBA3  +  e3[(DBB4)  cos  -  (DBB3)  Y^  cos  +  DBC3  +  DBD3 

ADS  =  DBA4  -  e,{DBB4)Y^  cos  +  DBC4  4  DBD4 

AE]  =  sin  e2(DBFl)  4  sin  eg(DBHl) 

AE2  =  e2(sin  <;6jDBF2)  4  Z^.  sin  0^  (DBFl)]  -  DBGl 

-  CjLsin  <A^  (DBH2)  -  Y^  sin  6^  (DBHl)]  -  DBll 

AE3  =  e2[sin  (DBF3)  4  Z^,  sin  (DBF2}]  4  DBG2 

4  e.[sin  9i.^{DBH3)  -  Y^.  sin  (DBH2)]  4  DBI2 

AE4  =  e,[sin  (DBF4)  -  Z ^  sin  (DBF3)]  4  DBG3 

-  ej(sin  (DBH4)  4  Y^  sin  (DBH3)]  4  DB13 

AE5  =  -Z^.cj  sin  6^  DBF4  4  DBG4  4  Y^Cj  sin  DBH4  4  0614 

hi  =  -Ld  -  ^’d 

h2  =  ®13  +  eg^d  + 

^13  =  “  ''9®13Md  ~ 

®II®13*^d  +  ®8\  +  ^8®l4''*d 
*’14  =  ^8«]4^1a  -  ^9*=13'\ 

^21  = 

h2  =  L^Mq^^d  ^  -  !^/SLd  - 

^23  =  -V’g^d  +  ^aMd>^j8  +  L^V’q  '  ‘^/s'^d^s  " 

^24  =  ®8^’a^/3  - 

hi  =  -  ^’^^'x/lxx) 

h2  =  -  (Ixx'I*x>NjS  4 

hs  =  ^/3'h  "  ®]4^/8  +  «8«i  1*^/3  ~ 
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^4  =  "  ®11®13^;8 

1’41  =  -N^Md^Ixz-^IxJ  -  MdL^ 

^42  =  “^(3®ll^d  7 

^3  =  +  N/S^gMa  -  «xz/IxAN/3  +  ^g^d^jS  +  L^N/3^n 

-  eu^aLja  - 

>^44  =  egM^N^  - 

Cu  =  1  -  WxzIxzVa,,!,,)] 

Cl2  =  -14  -  ^’q  -  -7  -  <Ixz''lzz)-9  + 

Cjg  =  enM^  -  +  ege^  ^(1^^' [^^)  -  ejgejg(I^j,''I^j^) 

“  ^^*z/^xx^^q®12  ~  -g^^q^^xz'^^zz^  “  ®9®12  +  -ll^U  “  -S^lO 
^14  ^  -T-H^^'q  "  -8-11-12  +  -9-10-13  +  -9-12^'^q  ~  *7-]l-13  ~  -8-10-14 
^21  =  -  N^(lxz'Ixx^ 

^22  =  ~  +  (Ixz-Ixx^N^Mq  + 


Cgj  ^ 

+  euLySM,  -  - 

14^P  +  -S-ll'^/3  - 

N^M^eg  +  N^eg 

]  ~  -1I-13^J|3 

^24  = 

-14^/8 

+  -8-1 1^/9  - 

■  -9^>Mq  -  -11-13 

L/3 

*^31  = 

-N^- 

«xz'Izz)4 

*^32  = 

^>Mq 

-  N^e^  -  - 

•  (Ixz'Izz^L^  +  +  ej^L^ 

^33  "" 

--lO-ia’-^  +  + 

N^e?  -  e.MqN^  + 

-8*1 0^(3 

"  -12''^q^/3 

+  L^'i2  t  'lxz''l„»',L^ 

n 

II 

--]0-13^^  +  eg-lO^'/S 

C41  - 

DBBl 

=  Md^ii 

C,2  =  DBB2  =  -L^ejo  +  e,iN„  -  +  -O^d^xz^Izz^  +  -12Md(Ixz'''IxJ 

-'^d-lO^Ixz'Ixx^  +  -]lLd<lxz''lzz^  +  K  ■*  -14^’d  ~  -7^0 


49 


AEDC-TR-7S-138 


C^3  -  DBB3  -  -L^ej^e^Q  -  -  e^CjgMa 

+  «9«]0^’d  “ 

C44  =  DBB4  =  -e9e]2MQ  - 

C5.  = 

^52  =  -^/3®ll*^d  -  +  ®9^'^d^'j3  “  ^/3-'''a^^z/^xx^  + 

-  -  ®iAL/3  - 

*^53  =  “  ^^xz’^Ixx^''^a^'j8  +  ®9'^^d^j0  +  "  ®H^^a^|3 

-  -  ^/3®l/''d 

^54  =  •^9'V^/S  - 

^61  =  ^’Id^/^  - 

^6  2  =  I-/3^d®lO  +  +  «7^''d^^/9  +  ‘^''d^^  -  ^'/g^d^lO 

+  L^Ma(r,,/IJ  -  L^e,2M^  + 


^63  =  ^/3^10^‘d  “  ®7^<a’^'/3  +  ®7'^d'^^  ~  ~ 

^64  =  ~®12'\^y3  " 

DBAl  =  — TerapaCcj  cos  o^) 

DBA2  =  62  cos  0p[Tempa  ~  Tempb] 

DBAs  =  Eg  cos  ii|j[Tempb  Z^,  —  Tempc] 

DBA4  =  62  ^o  ^vx 
DBCl  ==  -cos<^^e5ei[L^M^  + 

DBC2  =  -cos-^^e^eitM^L^  +  +  MdN^dxz^Ixx> 

-  M^egN^  +  -  L^N^en  +  ‘\*-/3«lll 


DBC4  =  cos  ^oe.e,(- iM^Cj^L^  +  'Vla^9^/3^ 
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DBD!  =  cos 

DBD2  =  CQs4>^e^e^[\\^L^{l^^  \,)  -  ^  +  ’^’a^'/3 

-  -  CjoL^N^  + 

DBD3  =  co5  6^e.e3[-M^L^ej2  + 


DBD4  =  coS(^^  _  M^L^e , ,) 

DBF]  =  Lj0(l,/-I,,)  + 

DBF2  =  L^e,2  +  L/gM./Ixz  1,.;  -  ^7)  +  ■"'/3"t  +  '^^ 

DBI-3  =  L^e,oc,3  +  L^Cis^’q  -  ’-j3®l2  “  zz'>  + 

-  V;  - 

DBK4  =  l-^e,oe,3  *  L.^e^^M,,  +  N^e,M,j  - 

DBG]  = 

DBG-:  r=  -s>n  ^^eoi;M„L^(r^/I^;  -  '^'’a®7^’/S 

-(-  N^M„  -  Fa'^/3°]0  +  \*-/3‘'l0^ 

DBG3  =  sin  ci  “  ''a^/S 

+  1 0  "  ^d^/3®lol 

dbg;  =  sin  ri„e2<M^ej2L^  +  \e-N^) 


-<lx,  “  ^d 


DBin 
DRI12 

DBII3  =  -l.^c,oC,3  -  L-e,2M  -  \l^e,3e,  +  -  e,\l,N^  + 


=  '-‘a^\2  +  ^‘<i 


F^MJL,  L.)  +  e, 


■  q  ‘  ^  X  7  *  7  //  a  1 .1 


“  Vl2''3  ^  ^8^d®l0 


Dlill;  =  ^ 
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DBll  =  -sin^i^egejCL^Nj^  -  L^K^> 

DBI2  == 

-H^d  +  MqL^Nri) 

DBI3  =  -sin  ^^egejCM^L^ejg  +  L^Md®13  ” 

- 

DBI4  =  -sin  (^^e5ej(M^L^ej3  -  IVI„N^eg) 


Static  Aerodynamic  Coefficients;  Body  Axis 


Cjvj  =.  Normal  force/ qS 
Cy  =  Side  (orce/qS 

=  Axial  force/^ 


Cg  =  Rolling  moment/qSb 
=!  Pitching  moment/qSc 

Cjj  =  Yawing  moment/qSb 


Static  Stability  Derivatives:  Body  Axis,  per  radian 


Cjj  =  dCj^/dn 
a 

Cy^  =  <9Cy/diS 
Cg^  =.dCg/a^ 


C  =  dC  /da 
C  =  dC  /dB 

np  n 


Dynamic  Stabiiity  Derivatives:  Body  Axis 


Cy  =  (9Cy/5(pb/2VJ 
Cy^  =  dCy''<9(rb/2VJ 

Cgp  =  dCg/.5(pb/2V^> 
Cg  =  dCg''<?{qc/2V^) 
Cg^  =  dCg/d{i8b/2VJ 
Cg  T=  ^Cg/(?(dc/2V  ) 


Cg  =  i?Cg/d(rb/2VJ 

G„  =  dC^/d(qc/2VJ 
q 

=  aC^/(9(rb/2V  > 

mm  <» 

r 

C.  =  5C^/d(pb/2VJ 
P 

C  =  dC  /diah/2V  ) 

m  •  m  cw 

a 

c„  =  ac„/d{pb/2vj 

p 


.  aC/«rb/2VJ 

-  5C/«^b'2VJ 

C.^  -  i;C„/S(«/2VJ 

C,  .  <)C/S<qc/2V^) 
q 
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Dimensional  Stability  Derivatives:  Body  Axis 


L,  -  (pV.Sb2'4I,^)C,^ 

L,  -  (pV,Sb/2l.^)C,^ 

Lp  -  (/>'»Sb=  '4I„)Ct 
P 

L,  =  (pSb2/41 

L,  =  (pV^Sb2/4I,y)Cg^ 

*^4^  =  (pSb^ '41  COS  d  )Cn 

o 

Ld  =  ^oL-w 

Lo  =  V  L 

P  DC  V 

cos 

M^.  =  (pSc2  '41  cosa^lC^  . 

^  ^  a 

M,  -  lpV.Sp2/4l„)C„^ 

'V  -  (pV„Sc2.'4l„.)C„^ 

M,  .  <pV.Sp^'4l,,,)C„_ 

Md  -  l-'.M, 


N.  -  VV.Sb'a^^C.^ 

N;  -  (pSb2'4l.,)C,^ 

Np  -  ((.V.Sl2/4I„)C„ 

P 

N,  .  (pV.Sb2'4I„)C„^ 

I'w  -  TO»4<.>C 

a 

N,  .  (pV.Sb2/4I,.)C. 

q 

NyS  = 

''''v  =  (p''«S''2m)C 

=  (pV^Sb/4m)CY 

^  r 

Yp  =  (pV_^Sb'4m)CY 
-Z^  =  (-pV_^S^2m  cos  +  C^) 
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NOMENCLATURE 


b 

c 

eg 

g 

Ixx ,  y  >  ^ZE 

Ixz 

it 


m 

n 


P,  q.  r 


P,  q,  r 


q 

s 

u 


V. 

X,  Y,  Z, 


Wing  span,  ft 

Wing  chord,  ft 

Aircraft  center  of  gravity 

Acceleration  due  to  gravity,  ft/sec^ 

Moments  of  inertia,  body  axis,  slug-ft^ 

Product  of  inertia,  body  axis,  s!ug-ft^ 

Aircraft  horizontal  stabilizer  incidence,  deg 
Aircraft  mass,  slug 
Aircraft  load  factor 

Roll,  pitch,  and  yaw  rates  about  aircraft  X,  Y,  and  Z  body  axes 

Roll,  pitch,  and  yaw  accelerations  about  X,  Y,  and  Z  body  axes,  rad/sec 

Dynamic  pressure,  Ib/ft^ 

Wing  area,  ft^ 

Velocity  along  body  X-axis,  ft/sec 
Free-stream  velocity,  ft/ sec 

Orthogonal  body-axis  system  fixed  in  the  aircraft  with  the  X-axis  along 
the  longitudinal  centerline  of  the  body,  the  Y-axis  normal  to  the  plane  of 
symmetry,  and  the  Z-axis  in  the  plane  of  symmetry 


SUBSCRIPTS 


a  Aircraft  angle  of  attack,  deg 

^  Aircraft  angle  of  sideslip,  deg 

p  Aircraft  density,  slugs/ft^ 

0  Aircraft  angle  of  roll 
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